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ABSTRACT 

Due to the lack of long term pulsed emission in quiescence and the strong timing 
noise, it is impossible to directly measure the braking index n of a magnetar. Based 
on the estimated ages of their potentially associated supernova remnants (SNRs), we 
estimate the values of the mean braking indices of eight magnetars with SNRs, and 
find that they cluster in a range of 1 ^42. Five magnetars have smaller mean braking 
indices of 1 < n < 3, and we interpret them within a combination of magneto-dipole 
radiation and wind aided braking, while the larger mean braking indices of n > 3 for 
other three magnetars are attributed to the decay of external braking torque, which 
might be caused by magnetic field decay. We estimate the possible wind luminosities 
for the magnetars with 1 < n < 3, and the dipolar magnetic field decay rates for 
the magnetars with n > 3 within the updated magneto-thermal evolution models. 
Although the constrained range of the magnetars’ braking indices is tentative, due to 
the uncertainties in the SNR ages, which come from distance uncertainties and the 
unknown conditions of the expanding shells, our method provides an effective way to 
constrain the magnetars’ braking indices if the measurements of the SNRs’ ages are 
reliable, which can be improved by future observations. 
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1 INTRODUCTION 

Pulsars are commonly recognized as magnetized neutron 
stars (NSs), sometimes have been argued t o be quark stars 
llZheng et al ] |2006l : 13^120071 : ILai et al.ll2013l l. A pulsar’s sec¬ 
ular spin-down is mainly caused by its rotational energy 
losses due to electromagnetic radiation, particle winds, star- 
disk interaction, intense neutrino emission or gravitational 
radiation (e.g., Peng et al. 1982; Harding et al. 1999; Al- 
par & Baykal 2006; Chen & Li 2006). An important and 
measurable quantity closely related to a pulsar’s rotational 
evolution is the braking index n, defined by assuming that 
the star spins down in the light of a power law 

- Kor , ( 1 ) 

where fl is the angular velocity, fl is the derivative of Q, 
and A is a proportionality constant. When the magneto¬ 
dipole radiation (MDR) solely causes the pulsar spin-down, 
the braking index is predicted to be n = 3. Note that this 
constant value of n = 3 holds only if the torque is propor¬ 
tional to D®, which is the case for the magnetic torque with 
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K keeping constant over time (Blandford & Romani 1988). 
However, in general n can vary over time because variations 
in any of these quantities may result in departures from this 
predicted value of 3. 

The standard way to define the braking index is 


uv ^ PP 

~ ^ i>2 - p2 ’ 


( 2 ) 


where D the second derivative of u — D/27r is the spin 
frequency, and P — l/v the spin period (see Manchester & 
Taylor 1977, and references therein). Be cautious that the 
measured value of n strongly depends on a certain time- 
span, and any variation in the braking torque will give a 
time-varying n{t). The braking index also reveals some sim¬ 
ple physical expectation. For example, n = 5 may point 
to gravitational-wave-emissior0 as the dominant spin-down 
mechanism, while in the case of n = 1, particle-wind radia¬ 
tion may be the cause of the pulsar spin down. 


^ For a normal radio pulsar, gravitational wave radiation can be 
only efficient during the first minutes/hours of the life (Manch¬ 
ester & Taylor 1977). 
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Table 1. Braking indices of 11 young radio pulsars. 


Source Name n Reference 


B0531-|-21(Crab) 

2.51(1) 

J0537-6910 

-1.5(1) 

B0540-69 

2.140(9) 

B0833-45(Vela) 

1.4(2) 

J1119-6127 

2.91(5) 

B1509-58 

2.839(1) 

J1734-3333 

0.9(2) 

B1800-21* 

2 (1) 

B1823-13* 

2 (1) 

J1833-1034 

1.857(6) 

J1846-0258 

2.65(1) 


Lyne et al. 1993 
Middleditch et al. 2006 
Ferdman et al. 2015 
Lyne et al. 1996 
Weltevrede et al. 2011 
Livingstone et al. 2007 
Espinoza et al. 2011 
Espinoza 2013 
Espinoza 2013 
Roy et al. 2012 
Livingstone et al. 2007 


*: For pulsars B1800—21 and B1823—13, the values of n = 2 
quoted in Espinoza (2013) are preliminary values of braking in¬ 
dices, with error bars of the order of 1 (Espinoza, private com¬ 
munication). 


For the majority of pulsars, it is extraordinarily diffi¬ 
cult to measure their stable and sensible braking indices n, 
from pulsar timing-fitting by using Eq.(2). The main rea¬ 
son is that the second derivative of v is very small and its 
detection is affected by the rotational irregularities, such as 
glitches, an abrupt increases in y and \y\ followed by relax¬ 
ations (Yuan et al. 2010), and timing noise (Zhang & Xie 
2013; Kutukcu & Ankay 2014; Haskell & Melatos 2015). 

By analyzing the timing residuals of 366 non-recycled 
pulsars over the past 36 yrs, Hobbs et al. (2010) mea¬ 
sured their braking indices and showed that they rang from 
—287986 to -1-3624. Unfortunately, the observed u for the 
majority of these pulsars with t > lO'^ yrs are dominated 
by the amount of timing noise present in the residuals and 
the data span. As noted in Pons et al. (2012), any small 
relative variation of the torque will produce large variation 
in n for old pulsars. In general, timing noise is ubiquitous 
among pulsars. Only for young radio pulsars, which spin 
down quickly and have stable first and second spin-period 
derivatives, their stable and reliable braking indices have 
been measured (shown as in Table 1), though these sources 
experience few small or large glitches, the post- glitch re¬ 
coveries dominate their timing behavior (Hobbs et al. 2010). 
Ten sources have positive braking indices below the canoni¬ 
cal value of 3. 

The very small braking index of n = 0.9 ± 0.2 (for PSR 
J1734-3333) and n = -1.5(1) (for PSR J0537-6910) could 
be consistent with the reemergence of magnetic fields af¬ 
ter a deep submergence into the star’s crust (Muslimov & 
Page 1995; Ho 2011; Vigano & Pons 2012; Gourgouliatos 
& Gumming 2014). According to the magneto-thermal evo¬ 
lution model proposed by Pons et al. (2012), a reemer¬ 
gence or a temporal rise of the crustal field, due to Ohmic 
dissipation and Hall drift, could in principle provide any 
braking-index value smaller than 3. Brak ing indices reach¬ 
ing v alues of ±10"^~® were also reported (|Manchester et al.l 
[ 200 ^. These observations co uld be contaminated b y unset¬ 
tled glitches or timing noise dAlpar fc Bavkall l2006l) . Other 
interpretations for such large |n| have also been proposed 
llBarsukov fc Tsvganll2010l : iBirvukov et al.l[2012l : IPons et ^ 
l2012l:IZhang fc Xiell2()12l) ^ 


Anomalous X-ray pulsars (AXPs) and soft gamma re¬ 
peaters (SGRs) are usually regarded as magnetars (Dun¬ 
can & Thompson 1992; Thompson & Duncan 1996). They 
are young NSs endowed with large surface dipolar magnetic 
fields. The soft X-ray emission in quiescence, bursts and 
flares of magnetars is supposed to be fueled by the rear¬ 
rangement and decay of the internal magnetic fields (e.g., 
Gao et al. 2011, 2012, 20130. 

It is commonly believed that pulsars, including tradi¬ 
tional radio pulsars and magnetars, are produced in core¬ 
collapse supernova explosions. Supernova remnants (SNRs) 
are the expanding diffuse gaseous nebulas resulting from the 
explosions of massive stars. SNRs can provide direct insights 
into the NS progenitor models and the explosion mecha¬ 
nisms, and thus have been studied by many groups (e.g., 
Vink & Kuiper 2006; Tian et al. 2007; Martin et al. 2014; 
Zhu et al. 2014). 

To date, we have obtained available information of 28 
known magnetars: 14 SGRs (11 confirmed), and 14 AXPs 
(12 confirmed). Among these sources, only 4 magnetars 
(IE 1547.0-540 8, PSR J1622-4950, SGR J1745-2900 
l Yan et all [2011) and XTE 1810-197) were discovered to 
emit radio waves. With the exception of SGR 0418-1-5729, 
all known Galactic magnetars lie within 2° of the Galactic 
Plane, consistent with a population of young objects, see 
“The McGill Magnetar Catalog” for relevant information. 
An online version of the catalog is located at the web site 
http://www.physics.mc gill.ca/~pulsar/magnetar/main.html 
m ausen fc Kaspill2014l) . 

The study of the braking indices of magnetars can pro¬ 
vide a wealth of information about the magneto—thermal 
evolution and the spin-down evolution of these objects. How¬ 
ever, due to the lack of long term pulsed emission in quies¬ 
cence and the strong timing noise, we have not yet mea¬ 
sured the values of n for magnetars. Observations hint that 
some magnetars have associations with clusters of massive 
stars. Recently, Tendulkar et al. (2012) measured the lin¬ 
ear transverse velocities, then estimated the kinetic ages of 
SGR 1806—20 and SGR 1900-1-14 assuming that the magne¬ 
tars were born in star clusters. Using the following implicit 
equation 


n = 1 -k 


TP 




( 3 ) 


and assuming Pi/P 1, the authors estimated n to be 
for SGR 1806-20 and l.lOtg'.gy for SGR 1900-kl4. 
Here T denotes the kinematic age of the magnetar (the time 
taken to move from the cluster to present position) and Pi 
is the initial spin period at t = 0 (Tendulkar et al. 2012). 
Although the braking indices of magnetars have been in¬ 
vestigated by some authors (e.g., Magalhaes et al. 2012), 
the work of specifically constraining the braking indices of 
magnetars has not appeared so far. In this paper, we try 
to constrain the braking indices of magnetars from the ages 
of their SNRs. In Sec. 2, we describe the potential mag- 
netar/SNR associations. In Sec. 3, we present constrained 
values of n for eight magnetars with associated SNRs, and 


^ The internal magnetic field of a magnetar could be a mixture of 
toroidal and poloidal magentic fields, and the toroidal component 
dominates (Mereghetti 2008; Olausen & Kaspi 2014). 
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possible explanations. In Sec. 4, we give a summary and 
discussion. 


2 POTENTIAL MAGNETAR/SNR 

ASSOCIATIONS 

It was predicted that about 10% of Galactic supernova ex - 
plosions may lead to magnetars jKouveliotou et al]ll994ll . 
Observations hint that about one third of the magnetars are 
likely to be associated with SNRs, suggestive of an origin of 
massive star explosions. Currently, the remnants ages are 
estimated based on the measurements of the shock veloci¬ 
ties, and/or the X-ray temperatures and/or other quantities. 
The SNR evolution can be schematically divided into three 
successive phases: 

(i) An ejecta-dominated phase. A shock wave created 
by the supernova explosion has a free expansion velocity of 
5000 - 10000 km s"^(Woltjer 1972). 

(ii) A Sedov-Taylor phase. The radiative losses from 
a SNR in this adiabatic expansion stage are dynamically 
Insignificant and can be neglected. 

(hi) A radiation-dominated snowplow phase. When the 
temperature drops below about 10® K, radiation becomes 
important, and energy is lost mainly through radiation. This 
phase ends with the dispersion of the envelope as its velocity 
falls to about 9km s“^ (Bandiera 1984). 

Sedov (1959) showed that there is a self-similar solu¬ 
tion for the adiabatic expansion phase (i.e., the structure 
of a SNR keeps constant within this stage), and firstly ap¬ 
plied this solution to estimate the age of SNRs in this phase, 
i.e., the Sedov age. The Sedov age estimates rely on many 
assumptions, but ultimately on the SNR size and the SNR 
temperature. If a SNR is too young (when the expansion 
is free), or too old (when radiation energy losses brake the 
expansion), the Sedov expansion is no longer applicable. 

Many authors (e.g., ©striker & McKee 1988; Truelove 
& McKee 1999; Gaensler & Slane 2006) subsequently inves¬ 
tigated the Sedov ages of SNRs. Based on XMM — Newton 
observations, Vink & Kuiper (2006) employed the Sedov evo¬ 
lution, and estimated the supernova explosion energies E, 
and the SNR ages fsNR by 

2.026i?tgNR = RIpo, (4) 

where po is the density of the interstellar medium (ISM), 
and Ra is the shock radius, which approximates to the SNR 
radius Rsnr. The shock velocity Va can be obtained by 

kT = ^ ^ 0.63mpUs, (5) 

where T is the post-shock plasma temperature, m-p is the 
proton mass, and the number 0.63 gives the mean particle 
mass In units of nip for a fully ionized plasma. From Eq.(5), 
we can obtain the value of Va, given the measured tempera¬ 
ture T. The age of the remnant can be estimated from the 
shock velocity similarity solution 

^ 2Rsnr 

Us = itSNR = -p--. (6) 

5rsNR 

Inserting Eq.(5) into Eq.(6), we get a convenient expression 
tsNR ~ 435 X RsnrT yrs , (7) 


where Rsnr and T are in units of pc and keV, respectively^. 
With respect to the potential magnetar/SNR associations, 
Gaensler (2004) proposed important criteria, which include 
evidence for interaction, consistent distances/ages, inferred 
transverse velocity, proper motion and the probability of 
random alignment. However, due to the observational limit, 
these criteria cannot always be applicable. The judgement 
as to whether a potential magnetar/SNR association meets 
these requirements ultimately depends on observational ev¬ 
idences. As pointed out by Gaensler (2004), one basic indi¬ 
cation of the likelihood of an association is the chance coin¬ 
cidence probability between a magnetar and its host SNR. 

In the following, we briefly review the measured and de¬ 
rived parameters for the eleven magnetars and their asso¬ 
ciated SNRs. All these SNRs have kinetic distances derived 
from their associations with other objects (e.g., molecular 
clouds) having known distances. 

(i) IE 1841—045 (hereafer IE 1841) The error box of 
this AXP lies within the SNR Kes 73 (G27.4-I-0.0), which is 
nearly circular in X-rays with an angular diameter of 0 = 
4.5 (Helfand et al. 1994; Gotthelf & Vasisht 1997). The 
probability of chance alignment was estimated to be 1 x 
10~^ for this association (Vasisht & Gotthelf 1997; Gaensler 
2004). Using the galactic circular rotation model, Tian & 
Leahy (2008) gave a range of revised distance to Kes 73 
about 7.5— 9.8 kpc, which implies a shock radius Rsnr ~ 
4.9 — 6.4 pc. Then the Sedov age of Kes 73 was estimated 
to be 500 — 1000 yrs for standard explosion energy {E = 
10®^ erg) and ISM n=lcm“® (Tian & Leahy 2008). 

(ii) SGR 0526—66 (hereafer SGR 0526) This SGR could 
be physically associated with the remnant N49 (Evans et al. 
1980) located in the Large Magellanic Gloud (LMG) at a 
distance of 50 kpc and interacting with a molecular cloud. 
The deep Chandra observation indicated the gas tempera¬ 
ture in the shell region to be r=0.57keV for N49 (Park et 
al. 2012). Assume an angular distance of 35 between the 
Shell region and the geometric center of N49, the SNR ra¬ 
dius is estimated to be 8.5 pc, then the Sedov age of the 
SNR is 4800 yrs (Park et al. 2012). However, the association 
between this SGR and N49 is controversial (e.g., Gaensler 
et al. 2001; Klose et al. 2004; Badenes et al. 2009), and more 
evidence is needed to testify this potential association. 

(hi) SGR 1627—41 (hereafer SGR 1627) This SGR is 
associated with the G337.0—0.1 (Hurley et al. 1999), and is 
located on the edge of massive molecular clouds. The angular 
size (radius of ~ 45 ) of G337.0—0.1 implies a radius of 2.4 
pc at a distance of 11.0 kpc (Sarma et al. 1997). The spurious 
probability of the SGR/SNR association was estimated to 
~ 5% (Smith et al. 1999). The age of the SNR was estimated 
to be 5000 yrs based on the Swedish-ESO Submillimeter 
Telescope observations (Gorbel et al. 1999). 

(iv) SGR 0501+4516 (hereafer SGR 0501) This SGR 
could be associated with SNR HB9 (also known as 
G160.9+2.6)(Gaensler & Ghatterjee 2008) with an angu¬ 
lar diameter of 128 (Leahy & Aschenbach 1995), but this 
association is not confirmed by observations (Barthelmy 
et al. 2008). The centrally brighten X-ray emission with 


® The propagated error of Rsnr is given by = 435 X 

/A^RsNH I ^l/2 

I- T -'-4t3- > 7“^®' 
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T ~ 0.82 keV from HB9 can be explained by self-similar evo¬ 
lution in a cloudy interstellar medium proposed by White & 
Long (1991) (WL91). Applying the WL91 evaporative cloud 
model, Leahy & Aschenbach (1995) (LA95) gave a distance 
of d =1.5 kpc (corresponding to i?sNR ~ 28 pc) and an age 
of 7700 yrs for the SNR. Based on the H I absorption mea¬ 
surement, Leahy & Tian (2007) (LA07) obtained a distance 
of d = 0.8±0.4 kpc (the mean radius is 15 pc at d = 0.8 kpc), 
and an age of 4000—7000 yrs for the SNR within WL91 evap¬ 
orative cloud model. In the following, we simultaneously con¬ 
sider both the results of LA95 and LT07, i.e., the distance is 
assumed to be 0.4 — 1.5 kpc, and the age for the SNR varies 
between 4000 and 7700 yrs. 

(v) PSR J1622—4950 (hereafer PSR J1622) The po¬ 
sitional coincidence of this AXP with the center of 
G333.9-1-0.1 with an angular diameter about 12 suggests a 
possible magnetar/SNR association (Anderson et al. 2012). 
G333.9-1-0.1 has a radius of 18 pc at the distance of 9 kpc, 
and is in a low density medium, suggesting the Sedov-Taylor 
expansion. (Anderson et al. 2012). Since no X-ray emission 
was detected from the SNR in the XMM-Neuton observa¬ 
tion (Anderson et al. 2012), the upper limit on the count-rate 
was estimated as 0.04 counts s“^ utilizing standard errors 
and roughly accounting (Romer et al. 2001). Assuming a 
thin-shell morphology and the Sedov-Taylor solution (Sedov 
1946, 1959; Taylor 1950), the upper limit on the count-rate 
implies an upper limit on the preshock ambient density to 
be no=0.05cm“® for an explosion energy of 10®^ erg. Then 
the upper-limit for the SNR age was estimated to be 6 kyrs 
(Anderson et al. 2012). 

(vi) IE 2259-1-586 (hereafer IE 2259) This AXP lies 
within 4 of the geometric center of CTB 109 (Fahlman & 
Gregory 1981) with an angular radius of 18. 5 ± 1. 0 as es¬ 
timated from the XMM-Neuton EPIC images by Sasaki 
et al (2004). The probability of chance alignment is about 
10“^ for this association (Gaensler 2004). There are some 
more recent estimates of the distance to CTB 109 and the 
AXP. Kothes & Foster (2012) collected all the observational 
limits of the distance (Kothes et al. 2002; Durant & van 
Kerkwijk 2006; Tian et al. 2010), and estimated a consensus 
distance of 3.2 ± 0.2 kpc to CTB 109, which places the SNR 
inside the Perseus spiral arm of the Milky Way. Castro et al. 
(2012) reported the detection of a GeV source with Fermi 
at the position of CTB 109. By simulating in the CR(cosmic 
rays)-Hydro-NEI (nonequilibrium ionization) model of Elli¬ 
son et al.(2007), Castro et al. (2012) derived an age of the 
SNR to be 8.5 — 15 kyrs, which is superior to the previous 
age estimates. Based on the Chandra data and the Sedov- 
Taylor solution, Sasaki et al. (2013) estimated an age of 
10 — 20 kyr for the SNR from a two-component model for 
the X-ray spectrum. Sasaki et al. (2013) also concluded that 
their (new) Chandra study is consistent with the Fermi 
results by Castro et al. (2012). 

(vii) CXOUJ171405.07-381031 (hereafer CXOU 

J1714) Halpern & Gotthelf (2010a) firstly suggested 
an association of this AXP with CTB 37B having a size of 
17 (see the SNR Catalog]), and a distance of 10.2±3.5 kpc, 
Caswell et al. 1975). This association is possible because 
the first TeV source HESS J1713—381 is hosted in CTB 

^ at http://www.physics.umanitoba.ca/snr/SNRcat/ 


37B, and is coincident with the AXP CXOU J1714 (e.g., 
Sato et al. 2010; Halpern & Gotthelf 2010b). Horvath & 
Allen (2011) gave an extensive review of this association 
and the SNR age. Clark & Stephenson (1975) obtained age 
estimates for both CTB 37A and CTB 37B of ~ 1500 yrs, 
and claimed that either source could be SN 393. Such an 
association is possible, but the SNRs may be considerably 
older (Downes 1984). Assuming an electron-ion equilibrium, 
Aharonian et al. (2008) gave a Sedov age of ~ 4900 yrs for 
CTB 37B. The inferred age for the SNR may be decreased 
to ~ 2700 yrs due to a higher ion temperature and efficient 
cosmic ray acceleration (Aharonian et al. 2008). Note that 
the electron-ion equilibrium is usually not achieved in young 
SNRs. By introducing “magnetar-driven” injected energy, 
Horvath & Allen 2011 obtained an age of ~ 3200 yrs (see 
Fig.l in their work) for CTB 37B, but the explosion energy 
and the ISM density were always fixed to be 10®^ erg and 
lcm“®, respectively. Based on the Suzaku observations, 
Nakamua et al. (2009) obtained an ionization age of the 
thermal plasma associated with CTB 37B as 65013 qq° yrs, 
which is consistent with the tentative identification of the 
SNR with SN 393. 

(viii) Swift J1834.9—0846 (hereafer Swift J1834) This 
transient magnetar is located at the center of the radio SNR 
W41 and the TeV source HSS J1844—087 (Kargaltsev et 
al. 2012). Such a cental location of Swift J1834.9—0846 cer¬ 
tainly supports the association between this SGR with W41. 
H I observations from the VLA Galactic Plane Survey gave a 
close distance of 4±0.2 kpc, and an average radius ~ 19±1 pc 
for W41 (Tian et al. 2007). The authors estimated an age of 
about 60 kyr, using a Sedov model with E ~ 0.75 x 10®^ erg, 
and an age of about 200 kyrs from complete cooling expan¬ 
sion (Cox 1972). 

(ix) The AXP IE 1547—5408 lies within 30 of the center 
of G327.24—0.13 with a size of 5 in radio (see in the SNR 
Catalog). The physical association of theses two sources was 
firstly proposed by Gelfand & Gaensler (2007), and was es¬ 
tablished by other observations (e.g., Joseph et al. 2007; 
Vink & Bamba 2009), and the distance to the SNR is quite 
uncertain. Unfortunately, since there is no X-ray emission 
from G327.24—0.13, we cannot estimate the age of this SNR 
by its Sedov evolution, and the published age estimate for 
this SNR from other methods has not appeared so far. Thus, 
we cannot estimate the braking index of IE 1547—5408 An¬ 
other magnetar AX J1845.0—0258 lies with the 5 diameter 
SNR G29.6-1-0.1, and the probability of chance alignment 
between the source and the SNR is shown to be 1.6 x 10“® 
(Gaensler et al. 1999). However, due to a long-term burst be¬ 
haviors, there is a lack of detection of the period derivative or 
dipole magnetic field for AX J1845.0—0258, so the braking 
index of AX J1845.0—0258 also cannot be estimated. The 
newly identified magnetar SGR 1935-1-2154 was proposed to 
be coincident with the Galactic SNR G57.2-I-0.8 (Gaensler et 
al. 2014; Sun et al. 2011). To date, no reliable age or distance 
estimate is available for G57.2-I-0.8, and for the evaluation 
of the braking index of SGR 1935-1-2154. Thus, the three 
magnetars above will not be considered in the following. 

Now, we list the parameters for the eight magnetars 
and their SNRs. All the SNRs in Table 2 have published age 
estimates. Among these objects, a few SNRs have several 
possible values of tsNR, measured with different methods at 
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Table 2. Persistent parameters for eight magnetars and their SNRs. Column one to ten are source name, spin period, spin-period 
derivative, surface dipole magnetic field, characteristic age, possibly associated SNR, SNR distance, SNR radius, SNR age and reference, 
respectively. Data in Columns one to six are cited from Table 1 of McGill SGR/AXP Online Catalogue, up to Feb. 6, 2015. Data in 
Columns seven to ten are mainly cited from the SNR Catalog at http://www.physics.umanitoba.ca/snr/SNRcat/ and from the Coolers 
Catalog at http://www.neutonstarcooling.info/references.html. In case where the timing properties vary with either phase or time, or 
there exists multiple recently measured values, the entries are marked with an asterisk (*) and compiled separately in Table 3. 


Name 

P 

(s) 

o 

1 

Sd 

10^^ Gauss 

Tc 

kyrs 

SNR 

t^SNR 

kpc 

-Rsnr 

pc 

*SNR 

kyrs 

Ref. 

IE 1841 

11.788978(1) 

4.092(15) 

7.0 

4.6 

Kes73 

7.5-9.8 

4.9-6.4 

0.75+0.25 

[1-3] 

SGR 0526* 

8.0544(2) 

3.8(1) 

5.60(7) 

3.36(9) 

N49 

50 

8.5 

4.8 

[4-5] 

SGR 1627 

2.594578(6) 

1.9(4) 

2.25(9) 

2.2(5) 

G337.0-0.1 

11.0 

2.4 

5.0 

[6-8] 

SGR 0501 

5.76209695(1) 

0.594(2) 

1.9 

15 

HB9 

0.4-1.5 

7.5-28 

5.85-1.85 

[9-11] 

PSR J1622* 

4.3261(1) 

1.7(1) 

2.74(8) 

4.0(2) 

G333.9+0.0 

9 

18 

< 6.0 

[12] 

IE 2259 

6.9790427250(15) 

0.048369(6) 

0.59 

230 

CTB 109 

3.0-3.4+ 

15-17+ 

15+5 

[13-15] 

GXOU J1714* 

3.825352(4) 

6.40(5) 

5.0 

0.95(1) 

CTB 37B 

10.2 

25.2 

1.75+1.4 

[16-18] 

Swift J1834 

2.4823018(1) 

0.796(12) 

1.42(1) 

4.9 

W41 

3.8-4.2 

18-20 

130+70 

[19-20] 


Footnotesd In Sasaki et al. 2013, a very close value oi d = 2.9 — 3.5 kpc was adopted, because the distance to the density peak of H I associated with arm 
was permitted to vary between 2.9 kpc and 3.6 kpc, although the distance to the SNR is mainly similar at 3.2 kpc, according to Kothes & Foster 2012. 

Refs: [1] Helfand et al. 1994 [2] Gotthelf & Vasisht 1997 [3] Tian Leahy 2008 [4] Evans et al. 1980 [5] Park et al. 2012 [6] Sarma et al. 1997 [7] Hurley et al. 
1999 [8] Corbel et al. 1999 [[9] Gaensler Sz Chatterjee 2008 [10] Leahy & Aschenbach 1995 [11] Leahy Tian 2007 [12] Anderson et al. 2012 [13] Kothes et al. 
2002 [14] Kothes &; Foster 2012 [15] Sasaki et al. 2013 [16] Caswell et al. 1975 [17] Aharonian et al. 2008 [18] Nakamua et al. 2009 [19] Tian et al. 2007 [20] 
Kargaltsev et al. 2012 


Table 3. Variable or alternative values of persistent pa¬ 
rameters for three magnetars with SNRs. All data are cited 
from Table 2 of McGill SGR/AXP Online Gatalogue at 
http: //www.physics.mcgill.ca/~pulsar/magnetar/table2.html 
up to Jan 10, 2015. 


Name 

P 

(s) 

P 

10 "++ s/s 

Rd 

10+"+ Gauss 

Tc 

kyrs 

SGR 0526 

8.0470(2) 

6.5(5) 

7.3(3) 

2 .0(1) 

PSR Jf622 

4.3261(1) 

0.94-1.94 

2.5(4) 

5(2).5 

GXOU J1714 

3.823056(18) 

5.88(8) 

4.80(4) 

1.03(1) 

— 

3.824936(17) 

10.5(5) 

6.4(1) 

0.58(3) 


different times, we select the age tsNR that is the latest and 
most reliable (e.g., tsNR ~ (15 ± 5) kyr for CTB 109, as 
shown in Table 2). For convenience, we denote each SNR 
age (Isnr) in the same form of “main value -I- age error”. 
With the exception of SNR W41 (tsNR ~ 130 ± 70kyrs), 
all the ages of the magnetar SNRs are mainly distributed 
within a rather narrow range of 1 — 20 kyrs. 


3 CONSTRAINING MAGNETAR BRAKING 
INDICES 

3.1 Constrained values of n for magnetars 

Assuming a constat value of n = 3, the characteristic age of 
a NS is defined as Tc = P/2P. As we know, Tc is usually to 
be used to as an approximation to its real age, with which 
it coincides only if the current spin period was much larger 
than the initial value, the magnetic-field strength B, the 
moment of inertia, and the angle between the magnetic and 
rotational axes have been constant during the entire NS life. 


By integrating Eq.(l), we get the inferred age as 


t = 


{n-l)P 
R 


t = 2rc In 




(n 7^ 1), 


(n = 1) , 


( 8 ) 


given a constant braking index n, where t is the inferred age 
(e.g., Manchester et al. 1985; Livingstone et al. 2011). 

In order to constrain the values of n for magnetars, we 
make the following three assumptions: 


(i) Firstly, we assume that a NS’s braking index keeps 
constant since its birth. This is a strong necessary assump¬ 
tion when estimating the age of a very young NS. In deed, 
one cannot derive an analytical expression, as Eq. (8), to 
estimate the age of a NS, if we consider either K or n in Eq. 
(1) is a time-dependent quantity. The measured braking in¬ 
dex is obviously influenced by the variation in the braking 
torque during the observational intervals. Nevertheless, the 
derived braking index can provide useful constraints on the 
realistic one. 

(ii) Secondly, we assume that the age of a SNR is the age 
of a NS. Zhang & Xie (2012) suggested that the physically 
meaningful criterion to estimate the true age of a NS is the 
NS/SNR association. This suggestion is surely in the same 
way applicable to magnetars with hosted SNRs. 

(iii) Finally, we assume that the initial spin-period of 
a magnetar is far less than the current spin-period, i.e.. 
Pi P, because magnetars are supposed to be formed 
from rapidly rotating NSs (e.g. Duncan & Thompson 1992; 
Thompson & Duncan 1996). 

Base on the above assumptions, the inferred age in Eq. (8) 
can be estimated as 


P _ V 

(n - 1)P (n - l)i> 


( 9 ) 
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Table 4. Constrained values of n for the eight megnetars 
with SNRs. The alternative braking indices are marked with an 
asterisk(*), and calculated from the data in Table 3. 


Source n Timing Reference. 


Dib Sz. Kaspi 2014 
Tiengo et al. 2009 
Kulkani et al. 2003 


IE 1841 

13±4 

SGR 0526 

2.40±0.04 

— 

1.82±0.06* 

SGR 1627 

1.87±0.18 

SGR 0501 

6.3±1.7 

PSR J1622 

>2.35±0.08 

— 

>2.6±0.6* 

IE 2259 

32±10 

GXOU J1714 

2.1±0.9 

— 

2.2±0.9* 

— 

1.7±0.5* 

Swift J1834 

1.08±0.04 


Esposito et al. 2009a, b 
Gogu§ et al. 2010 
Levin et al. 2010 
Levin et al. 2010 
Dib & Kaspi 2014 
Sato et al. 2010 
Halpern & Gotthelf 2010b 
Halpern & Gotthelf 2010b 
Kargaltsev et al. 2012 


IE 1841 SGR0526 

I A- 

! ^q-cxo 
ps^j'ieft 

T # 

!sGR1627 


I □ — AXPs 
\ o — SGRs 


-2 -1 


Log,o( fsNRM- ) 


Figure 1. Logio(tsNR/7‘c) versus Logio-B. The dot-dashed line 
corresponds to £snr = Fq. Blue circles and red squares represent 
SGRs and AXPs, respectively. Data for empty circles and squares 
are from Table 3. For some sources (e.g., SGR 1627), the errors 
in data-points are smaller than the size of the symbols. 


for Pi ^ p. Inserting the relation t = £snr into Eq. (9) gives 

p P . , 

n^l -= 1 + ^- . (10) 

I'tsNR PtsNR 

It should be noted that t in Eq.(9) is the upper limit of 
the inferred age due to an approximation in Eq.(8), and n 
denotes the mean braking index of a magnetar since its for¬ 
mation due to the assumption of a constant braking index. 
In the following Sections, for the sake of convenience, the 
mean braking index will be called “the braking index” in 
short. 

We can calculate the values of n for eight magnetars 
with SNRs by combining the data in Tables 2 and 3 with 
Eq.(lO), which are presented in Table 4. Since the first term 
in the right side of Eq.(lO) is a constant, the error of n is 
solely determined by the variations of — (or of 
the second term in the right side of Eq. (10), and the stan¬ 
dard error An of n can be expressed as the form of a stan¬ 
dard propagated erroiQ. From Table 4, all the values of n 
deviate from 1 and 3. Five magnetars have smaller braking 
indices, 1 < n < 3, while the other three have larger braking 
indices n > 3, which imply that neither pure PWR model 
nor pure MDR model can account for the actual braking 
of magnetars. In order to investigate the relation of Tc and 
tsNR, we plot Logio(tsNR/'rc) versus LogioB for nine mag¬ 
netars in Fig. 1. 

In Fig. 1, magnetars with n > 3 are on the left of 
the dashed line whereas those with 1 < n < 3 are on 
the right. The large discrepancy between tsNR and Tc of 
a magnetar can be naturally explained as follows: Inserting 
P/p = —v/v = 2rc into Eq. (10) gives Tc/Isnr = (n — l)/2; 
if n > 3, then Tc > tsNR, and a magnetar looks “older” than 
its true age (e.g., IE 1841). In this case, the larger the brak¬ 
ing index is, the larger the disparity between Tc and tsNR, 
as illustrated by IE 2259. If 1 < n < 3, then Tc < tsNR, 
and a magnetar appear “younger” than it is, the smaller the 
braking index is, the larger the age-disparity. 


® There is a relation between An and 




3.2 Case of n > 3 


The influence of the magnetic held evolution on the braking 
indices of radio-pulsars has been well investigated in previ¬ 
ous studies (e.g., Geppert & Rheinhardt 2002; Aguilara et 
al. 2008; Pons & Geppert 2007; Pons et al. 2009; Vigan6& 
Pons 2012; Vigano et al. 2012). In this subsection, we focus 
on three high braking indices (n=32±10, 13±4 and 6.3T1.7 
for IE 2259, IE 1841 and SGR 0501, respectively). In order 
to explain why the three magnetars possess large braking 
indices, we consider the influences of crust magnetic held 
decay on a magnetar’s spin-down, and refer to the works of 
Pons et al (2009) and Vigano et al. (2013). 

According to Pons et al.(2012), the spin evolution of an 
isolated NS is dominated by the magnetic held anchored to 
the star’s solid crust, and a varying crustal magnetic held 
yields the braking index 


,,^3-4Li^=3-4:^^=3-4ll, 


Bi V 


Bp V 


tb 


( 11 ) 


where Bd is the surface dipole magnetic held. Bp = 2Bd is 
the surface magnetic held at the pole, and tb = Bp/Bp = 
Bd/Bd is the timescale on which the dipole component of 
the magnetic held evolves. From Eq.(ll), a decreasing sur¬ 
face dipole magnetic held, Bd, results in a decaying dipole 
braking torque, then the ehective braking index increases to 
a high value of n > 3. 

In order to reconcile the measured high braking in¬ 
dices for three magnetars with their (estimated) surface held 
decaying-rates, it is necessary to introduce the Hall induc¬ 
tion equation describing the evolution of crust magnetic held 
in the following form: 


dB 

dt 


= -V X 


X {e-'B) + 

47rcr ATverie 


[V X (e^'B)] X B 


( 12 ) 

where a is the electric conductivity parallel to the magnetic 
held, e‘' is the relativistic redshift correction, rie is the num¬ 
ber of electrons per unit volume, and e the electron charge. 
On the right-hand side, the hrst and second terms account 
for Ohmic dissipation and Hall effect (which redistributes 
the energy of the crustal magnetic held), respectively. As to 
the above differential equation, Aguilara et al. (2008) pre- 
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sented an analytic solution 

_ Bp _1 Bp 

dt "CQhm Bp THall 

where rohm is the Ohmic dissipation timescale, and rnaii is 
the Hall timescale corresponding to an initial surface dipole 
magnetic field at the pole, Bp. During a NS’s early magneto- 
thermal evolution stage, corresponding to the crust temper¬ 
ature T > 10® K, the average value of rohm in the crust 
is about 1 Myrs (Pons & Gepport 2007). Since both the 
magnetic field strength and the density vary over many or¬ 
ders of magnitude in NSs’ crusts, the average Hall timescale 
is roughly restricted as rnaii ~ Tohm/)! — 10)Hi3, where 
Hi 3 is the dipole magnetic field at the pole in units of 
10^® Gauss (Pons & Gepport 2007). We choose rnaii ~ 
5 X 10^ — 5 X 10® yrs, similar as in Vigano et al. (2013), 
and rohm ~ 10® yrs for the three magnetars. 

To estimate the values of Bd for the three magnetars, 
we refer to Vigano et al. (2013): 

(i) By introducing the state-of-the-art kinetic coefficients 
and considering the important effect of the Hall term, Vi¬ 
gano et al. (2013) presented the updated results of 2D simu¬ 
lations of the NS magneto-thermal evolution, and compared 
their results with a data sample of 40 sources including mag¬ 
netars IE 2259, IE 1841 and SGR 0501. It was found that 
by only changing the initial magnetic fields, masses and en¬ 
velope compositions of the NSs, the phenomenological diver¬ 
sity of magnetars, high-B radio-pulsars, and isolated nearby 
NSs can be well explained in their theoretical models. 

(ii) These sources have best X-ray spectra by Chandra or 
XMM-Newton, from which the luminosities and temper¬ 
atures can be obtained, their well-established associations, 
and known timing properties with the characteristic age and 
the alternative estimate for the age. Similar as in this work, 
Vigano et al. (2013) selected their SNRs’ ages as the ac¬ 
tual age estimates for the three magnetars (for IE 2259, 
IE 1841 and SGR 0501, logr*, (yr) ~ 4.0 - 4.3, 2.7 - 3.0, 
and 4, respectively, corresponding to Tk are 10—20, 0.5 —1.0, 
and 10 kyrs, respectively). For generality, the Hall evolution 
timescales for the crustal fields have been arbitrarily selected 
as 10®, lO'^, 10®, and 5.0 x 10® yrs, respectively (Vigano et 
al. 2013). 

(hi) Gomparing the cooling curves with iron envelopes for 
Bp = 10^® Gauss, IE 1841 and SGR 0501 could be born 
with even higher magnetic field of a few 10^® Gauss, which 
can provide a strong hard X-ray (20—200 keV) emission with 
a total luminosity ~ 10®® erg s“^ (Vigano et al. 2013). 

(iv) Gomparing the cooling curves for the group of NSs 
with Bd ~ 1 — 5 X 10^^ Gauss suggests that IE 2259 may 
be born with Bp ~ 3 x lO^'^ Gauss, needed to reconcile the 
observed timing properties and the persistent soft X-ray lu¬ 
minosity. 

According to the above arguments, we assume Bp = 
3.0 X 10^® Gauss for IE 1841 and SGR 0501, and Bp ~ 
3 X 10^^ Gauss for IE 2259. Inserting the values of Bp, Bp, 
Tohm, and THaii iuto Eq.(13), we can obtain the estimated 
values of the surface magnetic field decay-rates dBp/dt. By 
solving Eq.(13), we obtain the relation of dBp/dt and the 
field evolution timescale t for the three magnetars, as shown 
in Fig. 2. From Fig.2, dBp/dt decreases with increasing 
t. Due to the smallest value of dBp, the values of Bp/dt 



Evolution time t (yrs) 

Figure 2. The relation of dBp/dt and t for three magnetars with 
n > 3. The initial values of Bp are 1.4 X10*^®, 3.8 X10^"*, and 1.18x 
10^^ Gauss for IE 1841, SGR 0501, and IE 2259, respectively. 



Log,„ t 


Figure 3. The relation of n and t for three magnetars with n > 3. 

for IE 2259 are lower than those of the other two mag¬ 
netars. Assuming t — Tk (Vigano et al. 2013), the current 
values of dBp/dt of three magnetars are calculated as fol¬ 
lows: dBp/dt = —(1.32 X 10^^ — 6.55 x 10^^ Gaussyr"*^ for 
IE 1841, dBp/dt = —(5.19 x 10® Gaussyr“^ for SGR 0501, 
and dBp/dt = —(4.76 x 10® — 2.44 x 10®Gaussyr“^ for 
IE 2259. SGR 0501 has the largest value of dBp/dt be¬ 
cause it has the smallest age t = Tk = 0.5 — 1.0 kyrs. 
Thus, the values of timescale tb = Bp/Bp can also be 
estimated as follows: for IE 1841, tb ~ 1.07 — 2.14 kyrs, 
corresponding to t = Tk = 0.5 — 1.0 kyrs, for SGR 0501, 
Tb = 73.2 kyrs, corresponding to t = rt, = 10.0 kyrs, 
and for IE 2259, tb = 24.8 — 48.4 kyrs, corresponding to 
t = Tk = 10.0 — 20.0 kyrs. 

Assuming the surface dipolar field decay is the only fac¬ 
tor influencing magnetar braking indices, by combining the 
values of dBp/dt with Eq.(ll), we plot the diagrams of n ver¬ 
sus log t for three magnetars, shown as in Fig.3. In Fig.3, the 
solid-squares denote the values of n obtained from Eq.(lO) , 
and the empty-squares denote the values of n obtained from 
Eq.(ll) and the age-values presented in Vigano et al. (2013). 
The detailed calculations and comparisons are as follows: 

(i) For IE 1841, when t = Tk = 0.5 — 1.0 kyrs, which 
is consistent with the SNR age estimated by this work, its 
braking index n = 11.6 — 20.2, which is comparable to the 
measured value of n ~9—17 (see Table 4) in the magnitude. 

(ii) For IE 2259, when t = Tk — 10.0 — 20.0 kyrs, which is 
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also consistent with the SNR age estimated by this work, its 
braking index n=20—40, which approaches to the measured 
value of n ~ 22—42 (see Table 4) in the magnitude. 

(iii) For SGR 0501, when t = Tk = 10.0 kyrs, its braking 
index n =3.8, which is apparently smaller than the measured 
value of n ~ 4.6—8.0 (see Table 4). This discrepancy may 
be caused by: 1 ) the initial dipolar magnetic field estimated 
is too low; 2 ) the age estimate is too large (compared with 
t = tsNR = 4.0 — 7.7kyrs in this work), and 3)the decay of 
the magnetospheric braking torque, e.g., the surface dipolar 
field decay cannot be the only factor influencing magnetar 
braking indices. Here, we insert t = tsNR = 4.0 — 7.7 kyrs 
into Eq.(ll), then get n ~4.1—5.0 (as denoted in the empty- 
circle in Fig.3), which is slightly larger than the value of 
n =3.8 by Vigano et al.(2013). 


Although we have attributed the high-value braking in¬ 
dices of the three magnetars to the decay of the crustal sur¬ 
face magnetic field, there is the possibility of the decay of the 
magnetospheric braking torque, which can also result in the 
braking indices n > 3. Below we will discuss the spin-down 
evolution for the three sources. 

SGR 0501 was firstly discovered by the Swift 7 -ray ob¬ 
servatory (Barthelmy et al. 2008). From the RXTE/PCA, 
Swift/XRT, Chandra/ACIS-S, and XMM-Newton/PN ob¬ 
servations spanning a short range of MJD 54700—54950, 
its magnitude of v decr eased from —2.03(8) x 10~^^ to 
-1.75(6) X 10-^®Hzs"^ dGogfis et al.ll2008l. l2010lf . If this 
decrease in v is solely caused by the decay of Bp, the mean 
surface dipolar field decaying-rate of dBp/dt is —(4.6 ± 
0.8) X 10^® Gauss yr~^ (for a rotator in vacuum. Bp ~ 
6.4 X 10^^(—Gauss), corresponding to the braking 
index n ~ 63.1 — 88.3 obtained from Eq.(ll). 

Recently, Dib & Kaspi (2014) presented a summary 
of the long-term evolution of various properties of five 
AXPs including IE 1841 and IE 2259, regularly moni¬ 
tored with RXTE from 1996 to 2012. Influenced by fre¬ 
quent glitches and magnetar outbursts, IE 1841 and IE 
2259 exhibited large fluctuations in their spin-down pa¬ 
rameters V and v. However, both AXPs obviously exhib¬ 
ited net decrease in z>. For IE 2259, its v decreased from 
-1.020(17) X 10"^'‘ to -0.9674 8 x 10~^^Hzs~^ d uring an 
interval of MJD 50355— 54852 llDib &: Kas^ l2014l l. which 


corresponds to n ~ 1404.2 — 2653.8 according to Eq.(ll). 
For IE 1841, its z> decreased from —2.9954(7) x 10“^'* 
to —2.944(11 1 X 10~^^Hz s~^ during a interval of MJD 
51225—55903( Dib fc Kas^l2014h . which corresponds to n ~ 


9.8-15.6. 


From the analysis above, apart from IE 1841, the es¬ 
timated values of n from the surface dipolar magnetic field 
decay are far larger than the values of the measured mean 
braking indices for magnetars IE 2259 and SGR 0501. We 
suppose that these large and rapid decreases in v were ex- 
plained as the effe cts of timing-noise torque and outbursts 
dGogfis et al.ll201Cl l. instead of as a secular dipole magnetic 
field decay. An actual decay of magnetospheric torque may 
explain the high braking index for a magnetar, but could 
result in an over-estimate of the surface dipolar field decay 
rate. 

Many authors (e.g., Harding et al 1999; Spitkovsky 
2006; Beskin & Nokhrina 2007; Tasng & Gourgouliatos 2013; 
Antonopoulou et al. 2015) investigated the external magne¬ 


tospheric mechanisms that can explain larger decrease in 
\v\ during the pulsar spin-down. The rotation of a pulsar 
induces an electric field accelerating charges off the star’s 
surface. The magnetosphere is expected to be filled with 
cascade plasma, which will screen the electric field along the 
magnetic field lines (Goldreich & Julian 1969). The mag¬ 
netospheric current, J, flowing in the magnetosphere and 
closing under the polar-cap surface, provides an additional 
braking torque, 

-F = \j[r,[l,Bp\]ds. (14) 

This braking torque on the star from the magnetospheric 
current is comparable to — Tmdr, the braking torque on the 
star for MDR model 0 From Eq. (14), a decreasing mag¬ 
netospheric current J produces a decaying —Tj, which will 
cause a decrease in v. For an net decay of toque due to the 
decrease in J, the relative decrease in the spin-down rate is 
estimated as ATj/Tj ~ AQ/Q = Avjv assuming the mo¬ 
mentum of inertial 7, and the angle between rotational and 
magnetic axes, x are constants. Note that, a precise calcu¬ 
lation of D (or v) depends on solving the Euler dynamics 
equations that are very complicate and uncertain (Beskin & 
Nokhrina 2007). 

It is also possible that a decrease in x can cause a de¬ 
crease in spin-down rate, as observed. During the long-term 
observations, these three magneta rs with SNRs exper i enced 
frequent outbursts and glitches. dGogfis et al. I l2008l. l 2010 l: 
iDib fc Kaspill2014ll . During post-outburst or glitch recover¬ 
ies, the stellar platelets move towards the rotational axes. 
The magnetic field lines anchored to the crusts follow the 
motion,and the inclination angle x decreases (Antonopoulou 
et al. 2015). A decrease in x will result in a net decay of the 
magnetospheric braking toque. Thus, the overall effect can 
be a decrease in v. Eor a small decrease of x, the relative de¬ 
crease in is estimated as ATdip/TMDR — Ax/tanx — Av/v 
in MDR model, or as ATfr/Tg ~ Axsin2x/(1 + sin^x) — 
Ar/r, here Tg is the braking torque for the force-free magne- 
tospher4I|( Spitkovsky 2006; Antonopoulou et al. 2015). If the 
external magnetospheric mechanisms (whether due to the 
change in the current or due to the change in the inclination 
angle) dominate the spin-down evolutions for three magne¬ 
tars, the order of magnitude of the relative decreases in the 
braking torque is about |AT/T| « \Ar/r\ ~ 10“^ — 10“^, 
as observed in IE 2259, IE 1841 and SGR 0501. 

In summary, we have quantitatively analyzed several 
physical processes (the decay of surface dipolar field, the 
decline of the current, and the decrease in the inclination an¬ 
gle), possibly causing n > 3 for magnetars. A common char¬ 
acteristic of these processes is that they can result in a decay 
of the braking torque on the star. Here we ignore the effects 


® In the simplest approximation for pulsar spin-down in MDR 

^2 q3 

model, the braking torque on the star is —Tmdr = —sin^x? 

here R is the stellar radius, and x is the inclination angle, i.e., 
the angle between rotational and magnetic axes. 

^ If the magnetosphere is filled with abundant plasma it can be 
considered in the force-free regime, where the induced electric 
field along field lines disappears everywhere except in the acceler¬ 
ation zones above polar caps and regions where plasma flow is re- 

g2 

quired, the braking torque on the star, —Tff = ^ 3 -(IH-sin^x) 

(Spitkovsky 2006). 
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of the stellar multipole magnetic moments on the magnetar 
braking indices, because the higher-order magnetic moments 
decay rapidly with the distance to stellar center. 


3.3 Case of 1 < n < 3 

In previous works, many authors proposed various mod¬ 
els to explain why the observed braking indice s of pulsars 
1 < n < 3, e.g., neutrino and photon radiation (IPeng et al.l 
Il982f l: the combination of dipole r adiation and the pro¬ 
peller torque app lied by debris-disk (lAloar fc Bavkaj 120061 : 
IChen fc LJ 120061 1 ; frequent glitches, as well as magneto¬ 
sphere currents (Beskin & Nokhrina 2007). In this work, 
the lower braking indices of 1 < n < 3 for five magne¬ 
tars with SNRs are attributed to the combination of MDR 
and wind-aided braking. Recently, the extended emission 
observed from IE 1547 and Swift J1834 are proposed to 
be magnetism-powered pulsar-wind-nebulas (PWNs), which 
may accelerate plasma to very high energy and radiate high- 
energy photons (see Kargaltsev et al. 2012, and references 
therein). The possible existence of magnetar PWNs lends 
a potential and direct support to our assumption of wind- 
aided braking. 

With the aid of strong magnetar wind, the closed 
magnetic-held lines are combed out at the extended Alfven 
radius ta, where the magnetic energy density equals the 
particle kinetic-energy density. The braking toque result¬ 
ing from tangential acceleration of outflowing relativistic 
plasma i s —Twind ~ MQ,r\ , w here M is the total mass- 
loss rate dManchester et al.lll985l l. The accelerated outflow¬ 
ing plasma will carry away the star’s angular momentum, 
and thus the spin-down rate of the magnetar will increase. 
The added braking toque due to a magnetar wind extends 
the effective magnetic braking “lever arm ”, causing the 
magnetar braking index to decrease, i.e., n < 3. 

In the scenario of wind-aided braking, the value of the 
effective braking index n of a magnetar depends on the 
competition of the magneto-dipole braking toque —Tdip and 
wind-braking toque —Twind- If iTvindl |2dip|, the value of 
n will be close to 1, and if iTdipl ^ | Twind |, the value of n 
will be close to 3. Unlike the persistent soft X-ray luminosity 
of a magnetar Lx, the persistent magnetar wind luminosity 
Lw cannot be obtained from observations. Moreover, it’s 
very difficult to determine Lw of a magnetar theoretically, 
due to the lack of a detailed magnetar-wind mechanism. 
The well-known magnetar wind-braking model proposed by 
Harding et al. (1999) implies a braking index of 


2 ;/ 

I V I 27’ 


(15) 


where we employ standard values of 7? = 10km and mo¬ 
ment of inertia 7 = 10“^® gcm^. Note that, by dehnition, the 
inferred value of n is always within the range of 1—3. In 
the above expression, n is an effective braking index, thus 
Lw is the effective wind-luminosity, and 73d is the effective 
surface dipolar Held. However, due to lack of the detailed 
information of the effective n, 73d and Lw, the importance 
of Eq.(15) has not been recognized in previous works. From 
the above expression, the magnetar wind luminosity can be 


estimated as 

6(3 - 


(16) 


By using Eq.(16), we estimated possible wind luminosity of 
magnetars, Lw, as listed in Table 5. Here we assume that 
the effective value of surface dipolar field to be the observed 
value of 73d for convenience. 

The rotation-energy loss rate is defined as Lrot = 
7f2H = Atv^Ii/O = —4:71^lPP~^, here 7 = 10’’®gcm^ for all 
sources. From Table 5, the value of Lx is obviously unrelated 
to that of Lrot, this is because magnetar X-ray emission is 
powered by magnetic field free energy. With the exception 
of SGR 0526, magnetars’ wind luminosities are bigger than 
their soft X-ray luminosities. In addition, IE 2259 possesses 
a relatively large ratio of Lx/Lrot ~ 303.57 (see in Table 
5), which is higher than those of all magnetars. The large 
ratio of Lx/Lrot could be served as a strong hint of surface 
dipolar field decay for this source, as discussed in Sec.3.2. 

Recently, Tong et al. (2013) proposed a multipole field- 
powered wind model for magnetars. The authors assumed 
that all the high luminosity magnetars with Lx > Lrot 
spin down via wind braking, and transient magnetars with 
Lx < Lrot spin down via magnetic dipole braking. Based 
on a simple assumption of Lw = Lp = 10®®ergs“^ (Lp 
is the escaping particle luminosity). They concluded that, 
for typical AXPs and SGRs, the dipole magnetic field in 
the case of wind braking is about ten-times lower than that 
of magnetic- dipole braking, and a strong dipole magnetic 
field may be no longer necessary for some sources and so 
on. However, when estimating 73d, one should be cautious 
taking into account effective braking indices. 


4 SUMMARY AND DISCUSSION 

Considering the possible SNR associations and their ages 
and assuming a constant braking index over the magnetar’s 
lives, we have estimated the average values of braking in¬ 
dices,for the eight magnetar, and discussed the possible ori¬ 
gins of the braking indices. Our main conclusions are as 
follows. 

(i) Based on possible timing solutions of magnetars and 
updated measurements of their associated SNRs, we esti¬ 
mate the braking indices n of eight magnetars, and find that 
they cluster in a range of n ~ 1 — 42. 

(ii) We attribute the higher braking indices, n > 3, of IE 
2259, IE 1841 and SGR 0501 to the surface dipolar field de¬ 
cay, or to the decay of the external magnetospheric braking 
torque. In the former case, the estimated dipolar field de¬ 
cay rates, based on the updated magneto-thermal evolution 
models, are compatible with the measured values of n for 
these three sources. 

(hi) We attribute the lower mean braking indices, 1 < 
n < 3, of other five magnetars (SGR 0526, SGR 1627, PSR 
J1622, CXOU J1714 and Swift J1734) to a combination of 
MDR and wind-aided braking, and estimate their possible 
wind luminosities Lw- 

At last, we remind that, in order to derive an analyt¬ 
ical expression for estimating the magnetar’s true age, the 
braking index is assumed to be a constant. However, the 
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Table 5. The relations among Lrot, a-nd of magnetars with different braking indices. Five sources in the upper part of Table 5 
have braking indices 1 < n < 3, while three sources in the lower part have braking indices n > 3. All the date Lx(2-10keV) in quiescence 
are cited from McGill SGR/AXP Online Gatalogue. 


Name 

-^rot 

(erg s-l) 

-^X 

(erg s-l) 

Lw 

(erg s-i) 

Z/W/Z/rot 

Tw/Tx 

Tx/Lrot 

SGR 1627* 

4.3(9) X 10®"^ 

3.6 X 1033 

(2.304±0.997)xl035 

5.36±2.58 

64.0±27.7 

0.08±0.02 

Swift J1834* 

2.1 X lO®"' 

< 0.84x1031 

(3.16±0.17)x1033 

15.03±0.78 

37583.3±1964.3 

<0.0004 

CXOU J1714 

4.51 X 10^4 

5.6 X 1034 

(1.53±1.52)x1033 

3.39±3.37 

2.73±2.71 

1.242 

-1 

4.15x10^4 

5.6x1034 

(1.15±1.57)x1033 

2.77±3.78 

2.05±2.80 

1.35 

-1 

7.41(35) X1034 

5.6x1034 

(5.36±4.26)x1035 

7.23±5.76 

9.57±7.61 

0.76±0.04 

SGR 0526 

2.87(7) xl033 

1.89 X 1033 

(4.13±0.55) xl033 

1.44±0.19 

0.002±0.003 

65.85±1.61 

-1 

4.9(4)X1033 

1.89x1033 

(2.75±0.53)x1034 

5.61±1.18 

0.146±0.028 

38.57±3.15 

PSR J1622 

8.3(5) xl033 

4.4 X 1032 

(1.47±0.40)x1034 

1.77±0.49 

33.41±9.09 

0.053±0.003 

-1 

(7.0±2.4)xl033 

4.4x1032 

(4.87±4.13)x1033 

0.69±0.58 

11.1±9.4 

0.063±0.022 

IE 2259-1-586 

5.6 X 1031 

1.7 X 1034 




303.57 

IE 1841-045 

9.9x1032 

1.84 X 1033 




185.85 

SGR 0501-1-4516 

1.2 X 1033 

8.1 X 1032 




0.67 


Footnote:* According to the “fundamental plane” (Rea et al. 2012) for radio magnetars with Lrot > Axi this source should have 
radio emission, however no radio emission was reported to date. 


inferred quantities (e.g., Bd and Tc) of a magnetar strongly 
depend on the current values of P and P (or of v and !>) 
that can vary with time. There always exist observational 
errors when estimating tsNR for magnetars. Thus, for some 
sources, their constrained values of n are subject to substan¬ 
tial uncertainties, and will be modified with future observa¬ 
tions. We expect that more observations of magnetars and 
their SNRs will help in further constraining the magnetar 
braking mechanisms. 
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